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Abstract. We obtain an estimate from below for the remainder in Weyl's law 
on negatively curved surfaces. In the constant curvature case, such a bound 
was proved independently by Hejhal and Randol in 1976 using the Selberg zeta 
function techniques. Our approach works in arbitrary negative curvature, and 
is based on wave trace asymptotics for long times, equidistribution of closed 
geodesies and small-scale microlocalization. 



1. Introduction and main results 

1.1. Weyl's law. Let X be a compact negatively curved surface of area A with the 
Riemannian metric {g%j}- We assume that the Gaussian curvature K(x) satisfies 

(1.1.1) - K\ < K(x) < -K\ 

at every point x € X . Let A be the Laplacian on X with the eigenvalues = Ao < 
Ai < A2 < . . . and the corresponding orthonormal basis {&} of eigenfunctions: 
A(pi — \i4>i. Let N(X) = #{\/Ai < A} be the eigenvalue counting function. The 
asymptotic behavior of N(X) is given by Weyl's law ( |Holj ): 

(1.1.2) JV(A) = ^-A 2 + i?(A), R(X) = 0(X). 

47T 

It shown in [Berj that on a nonpositively curved surface R(X) = O (jAr)< In the 
present paper we study lower bounds for R(X). As in [J-Pj . one of our tools is 
thermodynamic formalism for hyperbolic flows. 

1.2. Thermodynamic formalism. Let G* be the geodesic flow on the unit tan- 
gent bundle SX and let be the (one-dimensional) unstable subspace for G', 
£ G SX. The Sinai- Ruelle-Bowen potential is a Holder continuous function 7i : 
SX — > R which for any £ £ SX is defined by the formula (see [B-R] . [Sin] ) 

(1.2.1) no = j t 



lndet dG'Is- 

=0 



For any continuous function / : SX — > M one can define the topological pressure 
(1.2.2) P(f)=snp(h fi + J fdfi 
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where the supremum is taken over all G'-invariant measures /i and denotes the 
measure-theoretical entropy of the geodesic flow (see [Bow] ). In particular P(0) = h, 
where h is the topological entropy of the flow. It is well-known that for the Sinai- 
Ruelle-Bowen potential P(—Ti) — and the corresponding equilibrium measure 
(i.e. the measure on which the supremum is attained) is the Liouville measure /il 
on the unit tangent bundle: 

(1.2.3) V = / Hdfi L 

Jsx 

1.3. Main result. Recall that /i(A) = 0(/2(A)) for a function f\ and a positive 
function / 2 means that limsup A ^ oc |/i(A)|//2(A) > 0. 

Theorem 1.3.1. Let X be a compact surface of negative curvature. Then 

(1.3.2) R(X) = ((In A) P< ~" /2> V g > 0. 

As was shown in |J-P[ section 1.4], the power of the logarithm in (|1.3.2[) is always 
positive: 

r(-n/2) >J ^ 

h - 2K 1 

Moreover, if the curvature is constant, the bound (jl.3.2[) reads 

(1.3.3) R(X) = n ((lnA)5~ £ ) 

for any e > 0. The estimate (11.3.31) was obtained in |Ran| and, in a slightly stronger 
form, in |Hej| section 17], using the Selberg zeta function techniques. Our approach, 
based on the Duistermaat-Guillemin wave trace formula, thermodynamic formalism 
and semiclassical analysis, allows us to treat the variable curvature case as well. 

Theorem 11.3.11 agrees with a "folklore" conjecture that on a generic negatively 
curved surface 

(1.3.4) R(\) = 0(\ E ) Ve>0. 

Genericity is important, since on arithmetic surfaces corresponding to quaternionic 
lattices one can prove a much better lower bound i?(A) = (et) ( see |Hej | ) . 

Remark 1.3.5. One may compare Theorem 11.3.11 with the lower bound for the 
pointwise error term obtained in |J-P| : 

(1.3.6) R X (X)= l^( x )l 2 -^ =^(VA (lnA)^^ 1 - 5 ) Ve>0. 

Estimates (11.3.21) and (|1.3.6|) are independent, and, in particular, (|1.3.2p can not 
be deduced from ljl.3.6p . Indeed, cancellations may occur when i? K (A) is integrated 
over a negatively curved surface (for instance, according to the conjecture (|1.3.4p . 
\/A should cancel out in the generic case). Also, the sequence of A-s yielding the 
fi-bound (|1.3.6[) depends on the lengths of the geodesic loops at x (see |J-P1 section 
5.2]), and hence for each point x such a sequence is apriori different. 
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1.4. Wave trace asymptotics for long times. Consider the spectral distribu- 
tion 

oo 

(1.4.1) Tre(i) = ^cos(V%V) 

i=0 

which is the even part of the wave trace on X . 

To prove Theorem 11.3.11 we use a modification of the Duistermaat-Guillemin 
asymptotic formula for the wave trace [D-G] . Originally, this formula captures the 
contribution of a single closed geodesic, however for the proof of Theorem 11.3.11 
we need to take into account the contributions of all closed geodesies of length 
Tq < L < T(A), where To > is some constant (see Lemma f2 . 2 . 1 [) . and T(A) — * oo 
at an appropriate rate as the spectral parameter A — > oo. 

Let x(t, T) be a cut-off function 

(1-4.2) X {t,T) = {l-m)p{^ 

where p € <S(K) is an even, non-negative Schwartz function such that supp/3 C 
[-1, +1], and ip(t) € Cg°(K) with ip(t) = 1 when t G [-T ,T ] and ^{t) = when 
|*| > 2T . 

The long-time version of the Duistermaat-Guillemin trace formula is given by 

Theorem 1.4.3. Let T(A) — > oo as A — » oo with T(A) < elnA for e > small 
enough, and let L 7 (resp. L$ ) denote the length (resp. the primitive period) of a 
periodic geodesic 7. Then the asymptotics of the smoothed Fourier transform of the 
wave trace is given by 
(1.4.4) 

Tr e(t) X (t,T) cos Xtdt = V ^(XL,) X (L 7 ,T) 

^ eLsp X<T(A) Vldet(/-P 7 )| 

Here, V-y is the linearized Poincare map corresponding to 7. 

The proof of Theorem ll.4.3l is the most technically difficult part of the paper. We 
use semiclassical microlocal analysis on small scales, see section 3. In order to avoid 
the "accumulation" of singularities in the wave trace and to make the stationary 
phase method work, we separate the contributions of each closed geodesic using 
Lemma 12.2.11 Although the idea of working to (suitably scaled) Ehrenfest times 
T(A) ~ elnA is well-established (see |Ber| IZel2[ IFa]). rigorously separating out 
large exponential sums with T(A)-terms from the wave-trace using small-scale H- 
microlocalization appears to be a novel approach to estimating remainders in the 
negatively-curved case. 

The main term in the asymptotics proved in Theorem 11.4.31 is given by the sum 
of the principal wave invariants at each closed geodesic. As follows from Lemma 
12.1.61 to prove Theorem ll.3.1l it is sufficient to show that this sum grows at the rate 
given by (|1.3.2[) . First, we use the Parry-Pollicott equidistribution result |P-Pj to 
calculate the asymptotics of the sum not taking into account its oscillatory nature. 
Then we deal with the oscillations in the wave invariants (the difficulty is that 
oscillating terms may cancel out in the sum) using a "straightening the phases" 
argument based on the Dirichlet box principle (cf. |J-Pj ). 
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1.5. Lower bound for R(X) in higher dimensions. Lower bounds for the error 
term in Weyl's law on higher-dimensional manifolds are in fact much simpler than 
on surfaces. The reason is that the contribution from the "singularity at zero" to 
the remainder dominates the contribution from the periodic geodesies. 

We write /(A) ^> X k for a function /(A) if there exist a constant cq > and 
a number Ao, such that /(A) > coX k for any A > Ao- The following L 1 -estimate 
(which is stronger than an il-bound) holds for the remainder in Weyl's law on 
a rather general class of Riemannian manifolds containing manifolds of negative 
curvature. 

Theorem 1.5.1. Let X be a manifold of dimension n > 3, such that J r ^ 0, 
where r is the scalar curvature. Then, 

(1.5.2) ~ / \R(t)\dt^ A"" 2 , 

X Jo 

Theorem 1 1 . 5 . 1 1 is proved using the asymptotics of the Riesz means (see }Safj ) in 
section 15.11 One can also prove Theorem 11.5.11 using the standard t — > + heat 
trace asymptotics, see |J-P1 §2.1]. In order to get more refined information about 
i?(A) on negatively curved manifolds of higher dimension it is natural to study the 
oscillatory error term, i? osc (A), see section 15^21 

2. TWO AUXILIARY LEMMAS 

2.1. Smoothed Fourier transform of the wave trace. In the notations of 
section 11.41 let 



(2.1.1) fc(A, T) = l°° cos(Ai) Tr e(t)dt, 

J — oo ^ 

where, T — T(A) will be chosen appropriately later on. Substituting (|1.4.1[) into 
(|2.1.1[) we obtain 



(2.1.2) fe(A,T)=^ff A)T (v^), 

i=0 

where, for r > 0, 

(2.1.3) H XtT{ r) = r ^ cos(Ai) cos W = ^ X ~ r)) + P ^ + r » 



Here, 



T 2 

p(T(X-r)) 



*n J-oc 



is the Fourier transform of p and the 0(A~°°)-crror in (|2.1.3[) follows from the fact 
that p(X + r) = 0(X~°°) uniformly for r > 0, since p £ <S(K). Replacing the sum 
in (|2.1.2p by an integral, we get the following representation of fc(A,T): 

(2.1.4) k(X, T) = r P(T(X-r)) +P (T(X + r)) dN{r) = /« 
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Formula (|2.1.4j) plays a key role in our analysis. We shall also use the following 
notation: 

/■oo 

(2.1.5) / H x , T {r)dR{r) = k(X,T) 

Jo 

Note that the contribution of the main term in Weyl's law has been subtracted 
from k(X,T) to get k(X,T). 
We shall use the following: 

Lemma 2.1.6. Let R(X) = o((lnA) b ), b > 0. Then k(X,T) = o((lnA) b ) uniformly 
in T for T > To, where Tq is an arbitrary positive number. 

Proof. By the assumption of the lemma, for any e > 0, R(X) < e(l + lnA) b for 
large enough A. Consider the left hand side of (|2. 1.5ft : 

/>oo 

(2.1.7) / H x . T (r)dR(r). 

Jo 

Taking into account (|2. 1.3ft and integrating (|2. 1.7ft by parts we obtain 
(2.1.8) 

fT f°° fT f 00 

<\T)<— J \ p >(T{r-X))\(\n{l + r)) b dr + — J \p'(T(r + X))\(ln(l + r)) b dr. 

Since p' is Schwartz class, the second term of (|2. 1.8ft is O(l). Changing variables 
in the first term of (|2. 1.8ft , we obtain 



rp p<X3 

— J^\p'(Ts)\(ln(l + X + s)) b ds = 
e(lnA) b 




du < Ge(lnA) 6 



for some constant C > 0, where the last inequality again follows from the fact 
that p' is Schwartz class. Clearly, the constant C can be chosen uniformly in T for 
T > To > 0. Since e can be taken arbitrarily small, we get k(X, T) — o((ln A) b ), and 
this completes the proof of the lemma. □ 

Remark 2.1.9. Lemma \2. 1.61 is proved similarly to the results of |J-P|, section 2.2], 
see also jK]. One may also compare it to Sar, Proposition 3.1]. Sarnak's argument 
gives the lower bound R(X) — Q(\/A) for the Weyl error on a surface under the 
assumption that the geodesic flow G* has a fixed point set of dimension two for 
some t > 0. This condition holds, for example, when G* is completely integrable, 
but it is not satisfied on a negatively curved surface. 

2.2. Separation of periodic orbits. In this section we shall prove the following 
dynamical lemma: 

Lemma 2.2.1. Let X be a negatively curved surface and let f2(j,e) denote the 
e -neighborhood of a geodesic 7 in SX with respect to the Sasaki metric. Then there 
exist positive constants To, B and 5' (depending only on the injectivity radius inj(X) 
and the lower curvature bound K\ ) such that for any T > T the sets J? (7, e - BT ) 
are disjoint for all pairs of closed geodesies 7 on X with length T 7 G [T — 6' ,T]. 

Note that since there are exponentially many closed geodesies on X of length 
Lj 6 [T — 5', T], disjoint neighborhoods have to be of exponentially small size. 
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Remark 2.2.2. Here and further on we write J?(Y, d) for the ^-neighborhood of the 
set Y . It should not be confused with the O notation for the lower bounds. 

Proof. Choose B > 2K\, where K\ is the curvature bound and K\ is an upper 
bound for Lyapunov exponents, cf. (|4.1.4j) . Also, choose 5' < inj(X)/3, and let To 
be such that 2e~ KlT ° < 5' . Assume for contradiction that there exist two closed 
geodesies 71 and 72 with T — 5' < L~ ll < L l2 < T such that the corresponding 
neighborhoods intersect, and that the geodesies are not inverses of each other (note 
that due to the restriction L 7 E [T — 8' , T] in the conditions of the lemma, the 
geodesies cannot be integer multiples of each other unless they are inverses). 

Denote the geodesies on X by jj(t),0 < t < L 7 . , and their lifts to SX by 
7j(i) = (jj(t),jj{t j); Jj(t) C X is sometimes called a footprint of Jj(t). Without 
loss of generality we may assume that 

distsA- ( 72 (0), 71(C))) < 2e- 2KlT . 

Since for any < t < L l2 

dist S x(7 2 (i),7i(*)) = dist SJf (G t 72(0),G t 7l (0)) < 2e- 2K ^e K ^ < 2 e ~^ T , 
and hence 

(2.2.3) distx ( 72 (i), 71 W) < 2e- KlT . 

In other words, the entire geodesic 72 lies in the 2e~ /flT -neighborhood of the geo- 
desic 71 and vice versa. 

For convenience, we rcparamctrisc 71 and define 

/3i(s) := 7l (L lS /L 2 ), 0<s<L 2 . 
By triangle inequality and the definition of (3\ , 

dist( 72 (i),/3i(i)) <dist( 72 (i),7iW) + dist(7i(t),/3i(i)) < 



2e- Klt + t\l - j^j < 2e- KlT + 5' < 

Accordingly, for any < t < L2 there exists a unique shortest geodesic a*(s) 
in X connecting 72 (t) and Pi(t). We shall choose the parameter s € [0, 1] so that 
at(0) = Ts(t) andat(l)=/3i(t). 

Define the mapping s) : [0, L2] X [0, 1] — > X by the formula 

$(t,s) = a t (s). 

We claim that $ defines a homotopy between 72(t) and f3i(t). Indeed, Q(t, 0) = 
72 (t), 1) = f3i{t). Moreover, since both 72 and /3i have period L 2 , we have 

a (s) = Q i2 (s), V s e [0, 1], 

and so $(•, s) is a closed curve in X. Finally, $(£, s) is continuous since the function 
dist(72(£), (h(t)) is a continuous function of t. 

On the other hand, /3i is just a reparametrization of 71, hence 71 and 72 lie in the 
same free homotopy class, contradicting the fact that on a negatively-curved surface 
there is a unique closed geodesic in each free homotopy class. The contradiction 
completes the proof of the lemma. □ 
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Remark 2.2.5. Lemma [2. 2. II is proved for the tangent bundle. The tangent and the 
cotangent bundles can be identified using the Riemannian metric, and therefore 
Lemma 12.2.11 holds in the cotangent bundle as well. In this setting it will be used 
in the next section. 

Remark 2.2.6. The proof of Lemma l 2 .2.11 generalizes verbatim to higher dimensions. 

Remark 2.2.7. In the literature on wave invariants |D-G[ iDonl IZellj it is customary 
to choose a cut-off function in trace formulas in the length spectrum so that only a 
single length of a closed geodesic is contained in its support. Since we go to T(A) ~ 
In In A times, it is not enough to localize exclusively in the length spectrum. Indeed, 
when we localize around a geodesic(s) of length Li, error terms in expressions like 
(|3.2.2p are of the order 0(l/|£j+i — where ... < Li < < ... denote 

distinct lengths of closed geodesies on X (ignoring the multiplicity). Therefore, the 
error would be large in the presence of "near-multiplicities" . We can not control 
the gaps in the length spectrum on a generic negatively curved surface, and hence 
we have to localize in the phase space as it is done in the next section. 

3. Wave trace asymptotics and small-scale microlocalization 

In this section we give a proof of Theorem 11.4.31 which is quite technical. Let 
us note that assuming Proposition l3.2.H the subsequent sections can be read inde- 
pendently of section 3. 

3.1. Plan of the argument. We choose a parameter T = T(A) > satisfying: 

• (i) T(A) — > oo as A — > oo 

• (n) T(X) < em A. 

Here, e > is a small constant that will be chosen later on. For the applications in 
this paper, it will only be necessary to take T(A) of order In In A, so we will not be 
concerned with determining the best possible constant e > in (|3.ip . 

To localize the contribution to the wave trace and to k(A, T) from a given closed 
geodesic 7 of length L 7 < T it is important to microlocalize the wave trace Tr e(i) to 
a neighborhood of 7 and then sum over all the different 7's. The complication here 
is that we want to take into account the contributions of all closed geodesies 7 with 
To < L 1 < T(X). Since T(A) — > 00 as A — + 00, the number of these geodesies blows 
up. As a result, simply summing stationary phase expansions for each of the 7's 
(which is automatic when T = 0(1)) is impossible when one needs to work with such 
long period intervals. The way to deal with this is to microlocalize on neighborhoods 
of the 7's that shrink fast enough as A — * 00 (but not too fast) and then split up the 
time interval [To, T(A)] into short "windows" of fixed size 8' > 0. In this context, it is 
natural to work with semiclassical pseudodifferential and Fourier integral operators. 
The crucial dynamical result we need here is Lemma l2.2.1l Roughly speaking, this 
lemma says that there exist neighborhoods .J? (7, e~ BT ) of 7 of size e~ BT in phase 
space T* X with the property that, for appropriate geometric constant 8' > 0, 
no other periodic geodesic with period in the window [L 7 — <5',L 7 + 8'] intersects 
f2("f,e~ BT ). Since T(A) < em A, this clearly suggests microlocalizing the trace to 

e -BelnA _ ^- Be _ ne ig n b ornooc Js of 7. 

3.2. A reformulation of Theorem 11.4.31 For the purposes of the proof of The- 
orem ll.3.1[ it is convenient for us to reformulate Theorem 11.4.31 in the following 
equivalent form: 
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Proposition 3.2.1. Let T{X) — * oo as A — > oo wii/i T(A) < elnA. Then, for e > 
small enough, 

^ L» cos(Ai^) v(£~,T(A)) / 1 

(3.2.2) k(x,t(\))= y -i — + Of ^ 

Here, V*y is the linearized Poincare map corresponding to 7 and T) is a cut-off 
function defined by (|1 .4.2[) . 

The proof of Proposition 13 . 2 .11 is divided in several steps that are carried out in 
sections 3.4-3.8. 

Equivalence of Theorem 11.4.31 and Proposition 13.2.11 immediately follows from 
the following 

Lemma 3.2.3. Let T(X) -^00 as X — ► 00. TTien 

1 f 00 

(3.2.4) k(X,T) = —— Tr e(t)x(t,T)cos(Xt)dt + 



Proof. Combining (12. . (|2.1.4p . (12.1.5|) and using that p is a Schwartz function 
we get 

(3.2.5) k(X,T)= f ^Ip-cos{Xt)Ti e(t)dt- 

J — OO 

A f°° 

— J r(p(T(X - r)) + p(T(X + r)))dr = 
cos(At) Tr e(t)dt - — f°° rp(T(X - r))dr + 0((AT)-°°) - 

cos(Ai) Tr e(t)dt - ^Xp(0) + O (J^J + ©((AT)—). 

Therefore, one may rewrite k(A, T) as 
1 f 00 

(3.2.6) k(A, T) = - p(t/T)tp(t) cos(At) Tr e(*)dt+ 

J — OO 

- J_ p(t/T)(l - V(*)) cos(Ai) Tr e{t)dt - — Ap(0) + O f ^ 

Consider now the first term on the right-hand side of (|3.2.6|) : 
1 f 00 

(3.2.7) - / p(t/T)tp(t) cos(Xt) Tr e(i)di 
^ J —00 

It follows from the trace formula |D-G| that the contributions to (|3.2.7|) from the 
non-trivial periods t = L 1 7^ with < L 1 < To are 0(1/T). At the same time, 
the wave trace at t — has the singularity expansion [Zellj : Tre(t) = —-^t~ 2 + 
a\t~ x + a2 + . . . , where the leading coefficient can be computed by integrating the 
principal on-diagonal term of the parametrix for the wave kernel, see |J-P| section 
3-1]. 

Taking the contribution of singularity at zero into account, we obtain that (|3.2.7[) 
can be represented as 

(3-2.8) A. mX + of± 
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Therefore, the leading term in (l3~2~7j) cancels the A-term in (|3~2~6l) (cf. [EEl Le mma 
3.2.1] where a local analogue of such a result is established). To complete the proof 
of the lemma we note that by definition T) = p(t/T)(l — ip(t)). □ 

3.3. Preliminaries and notations. We now briefly recall the calculus of small- 
scale fi-pseudodifferential operators |D-S[ [Sj] that will be needed to carry out the 
various microlocalizations. In the following we use the notation H = A -1 . 

Given < 6 < |, we say that a(x,^;h) e C^(T*X x (0, Ho]) is in the symbol 
class S 7 s n {T*X) provided 

\d£dla(x,Z;h)\ < C a ,f3h- m - 5 ^+W\ 

Given a g S™(T*X), one can define the corresponding ft-pseudodifferential opera- 
tor Oph(a) invariantly in terms of the Schwartz kernel: 

(3.3.1) Op h {a){x,y) = {2-Kh)- n [ e—P^KA a (x, £; h) ((r 2 (x, y))d£. 

Here, exp : T X X — > X is the geodesic exponential map, £ G T*X, r(x, y) is ge- 
odesic distance between x and y and £ S Cq°(R) is supported in a ball B € (0) 
and equal to 1 in _B e / 2 (0) with e > sufficiently small (one can take here e < 

in](X,g). Such operators form a calculus with Op h {Sf)o Op h {Sf) C Op h (S™ +m '). 
Calderon-Vaillancourt L 2 -boundedness and the fi,-Egorov theorem also hold [Sj| 
section 2]. Moreover, let x,y € M" be the coordinates of the points 6 I in 
some local coordinate system (here we abuse notation slightly and denote points 
on the manifold and their coordinates by the same letters). By using the Tay- 
lor expansion — exp~ 1 (y) = x — y + 0(\x — y\ 2 ), to make the Kuranishi change 
of variables £ i— > (1 + 0(x — y))£ in (|3.3.1[) and integrating by parts in £, one 
can locally rewrite (|3.3.ip in the somewhat more familiar form Opf L {a){x^y) = 
(27T?i)- n J Rn e^ x -^/ h a{x,i\ H)(l + Oih 1 - 25 )) C(\x - y\ 2 )d£,. However, it is useful 
here to work with the invariantly defined operators in (|3.3.ip . and we will do so 
without further comment. 

We will also use the following geometric notations. Let M be the universal cover 
of X. The fundamental group of X is denoted by T. For u eT, let L u = L 1 be the 
length of the unique closed geodesic 7 on AC that corresponds to the conjugacy class 
[uj] C T. The lifted bicharacteristic curve in T*M projecting to 7 is denoted by 
7 C T*M. The standard cotangent projection map is denoted by tt : T*M — > M. 
We use analogous notation on T*X. We will also denote a periodic geodesic on X 
(resp. T*X) and a choice of lift on M (resp. T*M) generally by the same letter 
when the choice of lift is uniquely specified. Similarly, functions in C°°(T*X) will 
be identified with their lifts to C°°{T*M). 

3.4. Hadamard parametrix. One needs to get an asymptotic exponential sum 
formula for the leading term in (|3.2.4|) just like in the standard case where T(H) ~ 1. 
As we have already indicated, we will need to microlocalize the wave trace on small- 
scales with respect to h = A -1 , so it's useful to introduce the parameter h in the 
lifted wave operator. We do this by writing the real part of the wave operator on 

M in the form e(t) = cosi\/A = cost (^^^j ■ Here we denote by A the Laplacian 

on M and by e(t) the real part of the wave operator on M. 
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Even though introducing ft amounts to simply rescaling the Hadamard parametrix 
approximations E]y(t) to e(t) (see l3.4.1"l below). it's useful to think of Em(x, y, t; ft); 
as the Schwartz kernels of a family of ft- Fourier integral operators (albeit, with 

trivial dependence on ft) where h\f~K. is a classical ft-pseudodifferential operator 
[D-SL chapter 7]. Since inj(M) = oo, the Hadamard parametrix approximation 
Ejq{t, x, y; ft) is valid for allt>0 and is given by the following well-known formula: 
[Bir] : 

(3.4.1) E N (t,x,y;h)=h-i^- f°° e »W*«)-*]/* (f^ a^y)^ 3 h k ) dO, 

In the last identity (|3.4.1[) . d t is understood in the distributional sense. Also, for 
a E R, 0± denote the standard homogeneous distributions |Ho2[ section 3.2] and 
a k S C°°(M x M); k = 0, 1,2, .... It is also well-known (see |Ber[ section 39]) that 
for fixed y E M, U — Un € C N (M x M + ) and moreover, for any compact Y C M, 
there exists C\ , C% > such that 

(3.4.2) sup |e(t, x, y; ft) - E N (t, x, y;h)\< C 2 e Cl|t| h N . 

y£Y 

In (|3.4.2|) the bound is uniform for (x,y) E M x Y with d(x,y) < do, where do is 
any constant, and C s = C 3 (d , Y, N): j = 1, 2. Then [Ber llCdV] . one can take 

(3.4.3) E N (t,x,y;h) = ^ E N (t,x,uy;h), 

ugr ; L„<T(fi) 

to be the parametrix approximation to the wave operator e(t) = cosi\/A on X. 

Moreover, there exists an appropriate cutoff function r\ E Cq } (M;[0, 1]) (see 
}CdV| p. 94], |Ber|. Lemma 34]) with diam(supp rj) < 2diam(A") + 1 so that 

(3.4.4) Tre(t)= [ E N (t, x, x; ft) rj(x) dvol{x) + C(e c|t| H N ), 

JM 

for an appropriate C > 0. 
From now on we put 

(3.4.5) E N<u (t, x, y; ft) := E N (t, x, uy; ft), 
and so, 

(3.4.6) E N (t,x,y;H) := ^ E N ^(t,x,y;h). 

uer ; L„<T(s) 

Then, since \t\ < e\ lnft|, the RHS in (|3"X2"j) is 0(h N ^- Cie ^) and so, by choosing 
e < C-f 1 and N large enough, it follows from Lemma T3 . 2 . 31 that 

«(A,r) = -i- V (°° [ E N , LU (t,x,x;h) v (x) x (t;T(h))e^ h dvol(x)dt 

1 f°° C 

+ 12 E N ^(t,x,x;h) V (x) X (t;T(h))e^ h dvol(x)dt 

(3.4.7) +o{T(ny 1 ). 

Remark 3.4.8. Note that here and further on the sum over u> E T is finite due to the 
presence of the cut-off function %(i;T(ft)): the summation is taken over elements 
u ^ id that such that < T(ft). 
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The integral in (|3.4.7|) with the the e ±l '/ fi appearing has a total phase function: 
^Mjas.y) = {r 2 {x,y)-t 2 )9±t. 

Since 9 > 0, 

ft^i.^z.y) = -2t0±l, 

and for any a > 0, 

d? x (t;T(h)) = O a (i), 

after first integrating by parts in 9 for the k > 1 terms in (I3.4.ip . followed by 
repeated integration by parts in £ one gets that 

^ V /" f £^ w (i,a!,a!;n) X (t;T(h)) e- u / H dvol(x)dt 



2T(ft) 



\-i\ 



(3.4.9) +0(T(fi)- i ), 



where, u (t, x, y; ft) is defined to be 



(3.4.10) 

^(*)| fe"^)-''!/* ^a fc (.,^)^ fc - 1/2 ft fc ) ^(-2i0±l)^. 

From now on, we denote both the Schwartz kernel and the associated operator 
by E^ M . We note that the error term on the RHS of (|3.4.9p is the sum of the 
0(T(fi,) _1 )-error in (|3.4.7[) and an 0(ft°°)-piece that arises from inserting the cutoff 
function ip(—2t9 ±1). To see this, split the ^-integration into the regions where 
\9\ < 1 and 16*1 > 1. Integrating by parts in t over the first region contributes an 
O (ft 00 ) -error in (|3.4.9p . The same thing happens over the region where \9\ > 1 by 
noting that on this set 

| -2t9±l\- 1 < \9\- x 

and that T(ft) = 0(|lnft|). 

The variables (t, 9) are roughly-speaking dual and the presence of the cutoffs 
ip(—2t9 ± 1) will be exploited later on when we will need to interchange orders of 
integration in the expression for k(A,T). The analysis for each of the two integral 
sums in (|3.4.9| is the same and amounts to time-reversal. We will denote the first 
sum in (|3.4.9j) by k+(A,T) and the second by k~(X,T), 

The next step is to microlocalize on shrinking ft-scales near individual lifts of 
periodic geodesies. First, we microlocalize on shrinking scales near S*X. 

3.5. Small-scale microlocalization near S*X. We will need to exploit the sep- 
aration of geodesies in phase space given by the dynamical Lemma l2.2.1l Since the 
latter result applies to unit speed geodesies on S*X, we begin by ft-microlocalizing 
the trace k(A, T) near the co-sphere bundle S*X on small-scales ~ ft A ; < S < 1/2. 
The argument here is quite standard, but for completeness, we sketch the proof. 
Let d : T*X x T*X — > R be the distance function defined by: 

d T *x({xi,£i), (2:2,6)) = d T x((xi,€*), (%2,&)), 

where ^,£2 are dual vectors to £1,6 and drx is the Sasaki distance on TX. 
Lemma 12.2.11 extends in a straightforward way to the cotangent bundle endowed 
with such a distance function, see Remark [2.2.51 
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Let C G Co°(K) be a cutoff function with ((x) = 1 for x near and satisfying 
suppC C [-1, 1]. For < 5 < 1/2 we define 

x 5 (x : d;h) = t(h- 25 d 2 (( Xl s;),s*x)). 

Then, since p £ 5*(M), it follows that, modulo 0(A~°°)-errors, in k(X,T) one 
can sum over only the eigenvalues Xj satisfying l^ 7 ^" ~ A| < A e for any e > 0. 
Write h = A -1 and from now on, consider only such eigcnfunctions, <f>j]j — 1, 2, .... 
Since (hy/A - I) fa = Xj{h)fa; j = 1,2,3,..., where \\j(h) - 1| = O t {h 1 - t ) and 
the fi.-pseudodifferential operator P{K) = hyfK — 1 is fr-microlocally elliptic off 
S*X — {(x,£) £ T*X; |£| g = 1}, by a parametrix construction in the calculus 
Opn(Sg), it follows that 



(3.5.1) \\Op H (l- X 5 (x,Z;h))fa\\ L z = 0(h°°). 

Since £\ P( T W( X 3 + A)) = 0(X~°°), it then follows from (I2T41) and (13XT1) that 
«(A,T) = i^p(r(A)(A J -A))-^p(0) + O(T- 2 ) 

(3-5.2) = ~ p(T{X)(Xj A)) • (Op h ( X s )fa,fa) - ^p(0) + 0{T-% 

j 

In analogy with the construction of the small-cutoff function x S € Cq°(T*X) we 
choose C G Co°(^) as above with Q(u) = 1 near u — 0. We define 
X 5 (x,£;fr) = C(^ 2 ^ 2 ((a:^),^M)) • Xr{x), where Xi? € Cg°(M) with X fl = 1 
on supp(?7). Then, from Lemma r3.2.3l and (I3.5.2P it follows that: 

(3.5.3) x e ±ri ^x(^ dvol(x) dt + ©(T^)- 1 ). 

From now on, the cutoff x 5 wu l be included in all computations and for each 
u e T, the operators E^ u : Cg°(M) -> C£°(M) will be replaced by the microlo- 
calizations Op* h (x 5 ) ■ ■ Op h {x s ) ■ Cg°(M) -> Cg°(M). To simplify the writing, 
we will continue to denote the latter microlocalized operators simply by E^ . 

3.6. Small-scale microlocalization near periodic geodesies. The second part 
of the small-scale ft-microlocalization involves microlocalizing on shrinking ?i-scales 
near the lifts to T*M of individual periodic geodesies in T*X. In light of the 
previous section it suffices to take 7 C Q(S*M, H ) for any fixed S £ (0, 1/2). 

From Lemma 12.2.11 it follows that one can put disjoint tubular neighborhoods 
J7(7, e~ BT ^) around all lifts of periodic geodesies 7 with periods in the time win- 
dows [i, t + 6']; To < t < T(h). The constant B > is uniform and depends only on 
the curvature pinching condition (|l.l.ip . 

Fix the lift of a periodic geodesic 70 on M and choose e so that S := Be < h. 
Then, for (x, £) £ T*M, we define the small-scale cutoff functions 



(3.6.1) 



t* (x,t;h) = <;{h- 2S d 2 ((x,o, 7o ))) 
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Clearly, £ S${T*M), where we have arranged that 5 < 1/2. Then, for any 
u> £ r, we write 



E 



opU(L) ■ e% m ■ o Ph (d) + OpUcL) ■ e 



■ Ph (l 
Ph (l - 



C 



10' 

5 \ 



+o P * h (i - q )E% <u ■ o Ph (c° ) + o P * h (i - a /0 ) ■ E* tU ■ o Ph (i - c;, 

and taking traces of both sides, we split the RHS of (|3.4.7[) into three different 
integral sums. Taking into account that Tr A = TrA* and TrAB = TtBA, it 
follows that 

(3.6.2) 

where, 
(3.6.3) 

«n(A;T) 



: (A; T) = K ± (A; T) + 2k± (A; T) + n 22 {\; T) + 0(T(H)), 



(3.6.4) 



"12 



(A;T) 



2T(h) ^ 



E 



O P n(d ) 



OpUO]e ±lt/h x(t;T(h))dt, 



and 
(3.6.5) 



2T{K) 



2T(h) 



E 



Tr[E± ■ Op h (d )- PP»(1 " Ct)] e ±l * /ri x(i; T(ft)) dt, 



We now estimate each of the integral sums in (|3.6.3|) - (|3.6.5|) separately. Roughly 
speaking, one should think of the decomposition in (|3.6.3[) - (|3.6.5[) as follows: (|3.6.3|) 
gives the microlocal contribution of single periodic geodesic 70 to the trace, (|3.6.4|) 
consists of cross-terms which we will show are 0{h°°) and finally, (|3.6.5[) is estimated 
in the same way as (13.6.31) by successively microlocalizing around all other periodic 
geodesies 7 7^ 70 with L 7 < T(K). Due to the small-scale microlocalizations and 
ultimately, the splitting of the time scale into short time- windows, the expansions 
in k(X,T) arc no longer classical polyhomogeneous in h. For this reason, it is 
necessary to give a somewhat different argument than in the classical T(K) = 0(1) 
case (see |Donj ). 



3.7. Studying /^(A, T) in normal coordinates. The goal of this section is to 
prove the formula (|3.7.8[) . For (x, y) £ supp(?7) x supp(r;) C M x M, we have 

(3.7.1) E± t u ■ 0p h {O ■ °Ph&^ y> h ) = J^,<>' V, t; w, H). 
Here 

I±Jx,y,t;u>,h) = {2ttH)- 2 J J J e^- e ' 9+r ^ x ^ z ^- c ^ 1{ ^<^ h a N {x,z,e-u,h) 

x r)(x) C* (z, £ H) i>(-2te ± 1) ((r 2 (z, y)) dzd£d0, 

where apf(x, z, 8; u, h) = h~?0 + 2 J2k=o a k( x , uiz)h k 6^_ k . By symbolic calculus for 
S-pseudodifferential operators, £ S° satisfies ( |D-S[ p. 78] ) 

(3.7.2) Zl(z,S',h) = \&(z,Q\ 2 + h l - as &- 1 (z,bh) 
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with |C^ | 2 ,C7o _1 e Sg- Since inj(M, g) = oo, to simplify the writing somewhat, we 
put here (,{r 2 {z,y)) = 1 in the ft-pseudodifferential cutoffs (see (|3.3.f jl ). 

Since wcT acts by by isometries on (M, g), it follows that exp~l(ujy)-(duj~ 1 ) t t; — 
expj 1 (y) • £. and so the expression for I N u (x, y, t; u>, ft) above can be rewritten as 

(2irh)- 2 J J J e^^^^-^^yy^^^^ix^^^,^ 

x n(x) (z, £; ft) 4>{-2t0 ± 1) dzd£d0, 

Changing the variables uiz i— > z, (c&l> )*£ i— > £ and using that w € T acts by 
isometries, one gets 

4jx,y,t;u;,h) = (2nh)- 2 J J J e^^^^^^a^zAh) 

(3.7.3) x J7(a:) C* K) i>{-2tB ± 1) dzd£d6>. 

Here, w has been scaled out of the amplitude ajv so that ajy(x, z,9;H) — ft 2 9 + 2 ■ 

Let us now fix a global normal coordinate system centered at x £ M. With some 
abuse of notation we identify a point z £ M and a vector of its coordinates in this 
system: z = (^1,22) G T3M = R 2 . For instance, in the following z — toy means a 
vector (zi - (Loy)i,z 2 - (uy)^). 

Writing the Taylor expansion 

- expj^wy) ■ £ = (z - - (1 + 0(z - wy))£ 

and making the corresponding Kuranishi change of variables £ i— > £(1 + 0(z — toy)) 
in the above integral, we get a coordinate expression for (|3.7j) : 

I±JX, y, t; U>, ft) = (27Tfi)- 2 J J J jl-^+r^e^ym/H^ z & . ^ fi) 

(3.7.4) x ?7(x) (4 (a; -1 *, ft) V(-2i0 ± 1) dzd^6», 

where a' N {x, z, 9; u>, ft) = cln{x, z, 9] ft) • (1 + C(z — wj/)). Here u> (again, with some 
abuse of notation) is understood as the transformation R 2 — > R 2 , mapping the 
coordinates of z £ M to the coordinates of ujz £ M in the normal coordinate 
system centered at x £ M; du) is understood as the Jacobian of this mapping. Note 
that in formula (|3.7.5p and Lemma 13.7.61 below, u and du> are also understood in 
this sense. 

Recall from (13.7. 2p that one can write = |C* ? + h 1 ' 25 C^ 1 , where the second 
term Cw -1 e S° satisfies estimates of the form 

(3.7.5) 0?3£C&-V _1 *> = «,/3 (e cz "/r ff <H+IW) . 

Here, C denotes possibly different positive constants not depending on u>. The 
constants C may depend on a and /?, but this dependence can be ignored since we 
need to take into account only a finite number of derivatives: \a + 0\ < 2(m + 1), 
where m is given by (|3.8.4[) . 

To prove (|3.7.5p we first note that since C^, -1 G Sg, by the chain rule, differen- 
tiating the symbol gives the negative powers of ft. The estimate then follows from 
Lemma [3.7.61 below, since uj~ 1 z £ 7r(supp C^ -1 ) an< ^ therefore r(x,u~ 1 z) < C. 
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Lemma 3.7.6. Let M and uj be as above. Fix a normal coordinate system centered 
at x £ M, and let z € M be such that r(x,uj~ l z) < C. Let z — (z\,Z2), uj~~ 1 z — 
((w z)i, (uj^ 1 2)2) be, respectively, coordinates of z andui~ l z in this system, and let 
6 = (£ij&) G T z M and dujt t = ((dw'Oi, (^Oa) S T*_ lz M be the corresponding 
covectors. Then d%{uj- l z) = O a {e CL ") and d^{duj^) = O p (e CL "). 

Proof. Since uj is an isometry, we can identify the tangent spaces T X M = M. 2 and 
T UX M = M 2 using du>, so that for any point y £ M, the coordinates of the point 
y in a normal coordinate system centered at x coincide with the coordinates of ujy 
in the normal coordinate system centered at lux. Taking this identification into 
account we obtain 

(3.7.7) ((uj' 1 z) 1 , (u>~ 1 z) 2 ) =expJ*oexp a .(zi,z 2 ). 

As was shown in |Ber[ Appendix, Propositions 1 and 3] , the Jacobian of the ex- 
ponential map is bounded away from zero and the derivatives of the Jacobi fields 
have at most exponential growth at infinity. At the same time, the derivatives of 
the exponential map could be expressed in terms of the derivatives of the Jacobi 
fields jCh , p. 103], and hence also have at most exponential growth at infinity. 
Therefore, the derivatives of the map exp~* in ()3.7.7ll are 0(1) since by assump- 
tion of the lemma r(z,ujx) < C. The derivatives of the map exp^. in (|3.7.7p grow 
exponentially in r(x,z) < r(x,ujx) + r{iox,z) < r{x,tox) + C. At the same time, 
r(x,ojx) < L u + C by triangle inequality. This completes the proof of the first 
estimate in Lemma 13.7.61 

Consider now the ^-derivatives and let — (/3i , /3a)- Note that d^(duj t ^) = if 
Pi > 2 or /? 2 > 2. The first derivatives of cL/£ pull out components of du; 1 , which 
is the transpose matrix of the differential dco. Let us recall in what sense duo is 
understood here: it is the Jacobian of the map uj : M 2 — > M 2 given by (|3.7.7p . But 
we have already proved above that the derivatives of this map (in particular, the 
first derivatives that are components of the Jacobian) grow at most exponentially 
in L w . This completes the proof of Lemma 13.7.61 □ 

Let us go back to the formula (I3.7.4|) . By integration by parts in £ in (|3.7.4|) . 
modulo O(?i 00 )-errors, one can assume that r(z, coy) = r(oj~ 1 z, y) < ^. By carrying 
out the same argument as for the leading term |£* o | 2 below and using the derivative 
bounds (|3.7.5p . it follows that by taking 6 > sufficiently small, the contribution 
to k u (A,T) of the remainder term h 1 - 25 C,^ 1 (uj- 1 z , (duj)^; h) is 0{h c{ ^) for some 
C(S) > 0. From now on, with a slight abuse of notation, we ignore this remainder 
and rewrite |CL| 2 as Ct,- So, from (|3.7.4[> . we need to study the integral 



By integration by parts in the 9- variable in (|3.7.8p it suffices modulo 0(h°°)- 
errors to assume that for any fixed e > 0, 




' / / e^ t2e+r2 ^' z ^ z -^/ h a' N (x,z,9;uj,h) 
x r](x) Cl [uj^z, [dujft K) i){-2t6 ± 1) dzd£d6. 



(3.7.9) 



\r 2 (x,z) - t 2 \ < e. 



3.8. Stationary phase in (z, ^-variables and an expansion for «n(A,T). 
The goal of this section is to prove Lemma 13.8.91 We would like to apply sta- 
tionary phase with parameters in the (z, £)- variables in (|3.7.8p (see [Ho2] Theorem 
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7.7.5). For this, we need the following estimates for derivatives of the phase function 
4>{z,^;x, y,uj) ■= 9r 2 (x, z) + (z — ujy)^. First, we have that for (u>~ 1 z, g?o/£; x, y, uj) <E 
supp(C^ o ) x supp(r?) x supp(ry) x {lj;^ < T{h)}, 

(3.8.1) d?dl<l>(z,S;x,y,w) = O a ^{T{Kf). 

The estimate (P7HTT]) follows from the fact that T(h)- 1 < 9 < T^ 1 on supp ip(-2t9± 
1) and the following bounds: 

d^(z,^x,y,oj) = d2(6r 2 (x,z)) = 9d^r 2 (x,z) = 9d«(\z\ 2 ) 

= O a (9r 2 (x,z)) = O a (l)T(K) 2 ; |a|>0. 

Here, we have used (|3.7.9[) as well as the fact that r 2 (x, z) = z\ + z\ since the Zj's 
are geodesic normal coordinates at x € M. The mixed (z, ^-derivatives of </> are 
pointwise 0(1). Finally, 

d^(z, £ x, y, u) = z - cuy = 0(L U ) = 0(T(H)), 

where the last estimate follows from the triangle inequality for the distance function. 
Moreover, the O aj/ 3(l)-constants appearing on the RHS in (|3.8.1[) are all uniform 
in lo e r. 

The required lower bound on the norm of the gradient says that for (lo~ 1 z 1 duo 1 ^, x, y 
supp(2 o ) x supp(?7) x supp(?y) x {lo; L u < T(H)}, there exists another constant C > 
(uniform in all parameters including lo), such that 

(3.8.2) \d z <f>\ 2 (z,S;;x,y) + \d^\ 2 (z^;x,y) > + d z r 2 (x, z)9\ 2 + \z-ooy\ 2 ), 

Indeed, the lower bound in (|3.8.2p follows by Taylor expansion around the critical 
point £ = — d z r 2 (x, z)9, z — uoy and the Hessian lower bound {d 2 ^4>) ^> Id. The 
latter follows from the fact that d 2 z <j) = 9d 2 z r 2 (x,z)) = 9d 2 (z 2 + z\) = 0{9) = 
0(1), d z d ( 4> = 1 and finally, d%<j> = 0. Given (|gXT]) and ^EM, by Hormander's 
interpolation proof of stationary phase (|Ho2j Theorem 7.7.5), 

N 

e i[-t 2 +r 2 {x,uy)]0/H 
-1/2-fc 



I±Jx,y,t;h) = (2nh)- 1 / 2 J2 I 

k=o Jo 

x b k N (x, y, 9; lo, h) ^{-2t9 ± 1) 9+ 1/z ~ K d9 



Set 

(3.8.4) 

In (I3.8.3P we have 



1 - 25 



b k N (x,y,9;u;,H) = h k f] [a k (x, z) ■ (^(lo^z, dw^W 



£=—d z r 2 (x ,z)9 ,z—ujy 



(3.8.5) +0(h k+m - 1 T(h) 2 sup \d^(a k -CL)\). 



|o!+/3|<2(m+l) 



Let us now estimate the contribution to k^(A,T) coming from the error in 
t|3.8.5[> . Since it is well-known that d z ak = O a (e CLu: ) ( [Ber| Appendix, Lemma 4]), 
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we have that 

ft" 1 " 1 SUp \&?d?(a k • Ci)| « e C£« ft -2(m+l)«+m-lj 

|a+/3|<2(m+l) 

(3.8.6) < e Ce|lnR| ft (m+l)(l-2<5)-2 < fi (m+l)(l-2*)-2-Ce ^ tf-Ce 

Consequently, from (|3.8.3() . (|3.7.8| and (|3.6.3p it follows that the contribution to 
K^A, T) coming from the remainder in (]3.8.5[) for each k = 0, N, is bounded by 

(3.8.7) 

r ( fi ) 2 f°° f°° \M-o+a+i\\^r+.rrti\\\a. t ^it/Kiti2-Ce *k a -i/z-k 



T(h) 

T <L Ii ,<T(?i) 



/oo />oo 
/ |^(-2t6±i)|x(t; T{h))\d t e ±u ' h \h 2 - Cf - h k e+ 1/2 - k dddt 

,T(ft) ,1/To 

«ft fc+1 - C£ T(ft) Y, 0j 
t,..<t<k\ Jt ° Ji/n*) 



rT(ft) /-1/To 



L^<T(h) 
r T(h) 



f t- W 

^ h k+X-Ce T ^ e C\T(H)\ / |T(ft)|- 1 /2+fc dt « ftfc+l-Ce | T(fi) |3/2+fc_ 

In the second to last estimate in (|3.8.7p we have used the exponential bounds for 
the growth of the number of periodic geodesies |M-Sj . 

Choosing e > small enough and taking into account that T(K) < e\\nh\, it 
follows that the error terms for k > in (|3.8.5p are 0(h ai ) for some ol\ > 0. We 
are only interested here in working up to such errors. So, it is enough to consider 
only the principal term in (|3.8.5|) . Next, note that for each j > 1, 
(3.8.8) 

and so again after possibly shrinking e > in the T(h) < elnft, it suffices modulo 
0(h a2 ) to restrict the analysis to the j = case. By the same token, it suffices to 
restrict to the k = case in (|3.8.3p . 

To get rt 11 (A, T), we put y = x in (|3.8.3[) . integrate over x £ supp(?7) and sum 
over ui £ r;T < L u < T(H). Substituting the formula (|3.8.3|) in (|3.6.3|) and taking 
into account the estimate in (|3.8.8[) . for appropriate > 0, we get the following 

Lemma 3.8.9. 

4i(A,T) = -i- Y, e"I-* , +^(- w -)l/%(x,*,fl; W| fc)(W + )-i 

xr}(x)ih(-2t6±l)^- \e ±lt/h x(t;T(h))} dtdvol(x)d9 + 0(h a3 ) 
at L J 

/•OO /• />00 

Y e i9 [-* 2+r2 ^^) ±4 ]/%(x,x,0; W ,ft)(M + )-3 

, JO JM J-oo 



2hT(h) 



If) JM J-oo 

xr](x) i/)(-2t6 ± l)x(t; T(h)) dtdvol(x)dd + 0(h aa ). 



In the last line in (|3.8.9j) , we have used that the leading term comes from applying 
d t to the exponential e ±u l h since dtx(t\ T(H)) = d t [p(t/T) (1 - ip)(t)] = 0(1). 

For a fixed u> £ T, every term in (|3.8.9|) can be represented as a sum using 
the expansion (|3.8.5|) for b%. Let us estimate each term in this sum separately. 
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Given the estimate in (|3.8.8| the j = - term dominates in (13.8.5|) and so, up to 
0(h ai ) -errors, it suffices to estimate: 

(3.8.10) J J / ^[-f^M/^/^^ .^f . g^^^l^fl) 

xrj(x)a (x,ujx) {h9+)"^ x(i;T(7i)) ip(-2t6 ± 1) dt dvol{x) dd. 

3.9. Splitting into short time-windows. The goal of this section is to prove 
Lemma 13.9.61 We split up the iterated time integral in (|3.8. 10|) into a sum over 
'time- windows' of size j^y- The reason for this is that Lemma [2.2.11 only controls 

the splitting of geodesies in a time interval of size (note that the original time 
variable t has already been rescaled to j^f^ at this point). Consider a covering of 
[T , T(h)} by open intervals Ij of length 5', j = 0,... [ T( ^,~ r ° ], and let rjj € Cfi°(R) 
be a partition of unity subordinate to this covering. For j = 0, [ T 1, T ° ], we 
will need the following cutoff functions: 

Xj (t;T(h)) = r 1j (T(n)t)-x(t;T(h)) 

(3.9.1) =r l3 {T{h)t)-p{t/T{h))-{l-m)- 

Since we have already inserted the cutoff function ip(—2t6 ± 1) in (|3.8.10[) this 
allows us to apply Fubini and do the (t, ^-iterated integrals in (|3.8.10p first and 
the ^-integration last. Indeed, we just rewrite the total phase in (|3.8.10p as a sum: 

(3.9.2) r(x,wx) + ^ ± (t,6 
where, 

$ ± (i,6»;a;,w) = (±t - r{x,ux)) - 9 (t 2 - r 2 {x,ux)). 

At this point, we would like to do stationary phase in (t, 9), treating (x, w, K) as pa- 
rameters. Just as in (|3.7.8p we need to establish a couple of estimates for derivatives 
of the phase < i> ± . The first estimate 

(3.9.3) dfidl^^Oyx^) = O a>0 (r 2 (x,ujx)) 

with O aj/ 3-constants uniform in all parameters is immediate. Since dtde<& ± = 
— 2t, d 2 ^ — and d$ $ ± = —29 one gets the following lower bound for the (i, 9)- 
Hessian of $ ± : 

(dl e ^Y-(dl e ^)>^id, 

and so, by Taylor expansion around the critical points t = ±r(x,ujx) and 9 — 
— ri one gets the uniform lower bound 

(3.9.4) |9 t $±| 2 + \d e ^ | 2 > 1 (\t±r(x,u > x)\ 2 + \9 - 
So, by |Ho2] Theorem 7.7.5, it follows that for the expression in (|3.8.10p : 

x (h9 + )-i X (t] T{h)) ^{-2t6 ± 1) dt d9 dvol(x) 

= 2nh 1 / 2 [ e ir <»^/fya;)ao(:c,wz)-& (x,-du>* • d g r(x,z)\ g=ux ) ■ r(x,u>x)-i 
Jm 

(3.9.5) xi/>(±r(x,u)x)/T)(l-ip)(±r(x,ujx))dvol(x) + 0{h 2 ~ 2S e CL "). 
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In the last line of (13.9.5j) , we have again used the exponential bounds for the deriva- 
tives of ao [Berj as well as the uniform lower bounds for ^ in (|3.9.4| . The contri- 
bution of the 0(h 2 ~ 2S e CL ") error term to Kn(A,T) is then 

« E « ^- 25 e c l T ^)l|T(ft)|- 2 . 

^ ' u\L u <T{H) 

After possibly reducing the size of e > further, this term is then 0(h a5 ) for some 
a 5 > 0. 

We summarize what was shown so far in the following 



Lemma 3.9.6. 



, T(W)-Tq . 

11 [n > idjtu&r j=o Jm 

■ Xj(±r{x,ojx);T(h))r(x,u;x)-i T)(x) dvol(x) +0{h as ), 

for some a§ > 0, where r\ G C^°(M; [0,1]) is an appropriate cutoff function with 
diam(supp 77) < 2diam(X) + 1. 

3.10. Stationary phase in the x-variables. The last step involves expanding 
each term in the u-sum in (|3.9.6[) separately. For fixed id 7^ u) G T, one carries out 
stationary phase in (|3.9.6|) in the x- variables transverse to the lift of the periodic 
geodesic on X given by 

(3.10.1) 7 (w) supp(r?) C M;V x f u (x) - 0}, 

where fu(x) = r(x,u>x) is the displacement function. 

Since the dynamical Lemma 12.2.11 only controls separation of geodesies in time 
intervals of size jrcj) we estimate the summands: 

(3.10.2) ^ J2 [ e^'^Waofr^.^-drf.dirfaz)]^) 



2T(h) 



x Xj{^ E - r { x i ujkx )'iT{h))r{x,ujx) 2 r](x) dvol(x); j = 0, . 



, T(h)-To 
1 5' J 



separately. For fixed id 7^ u G T and j G {0, [ T ^r T ° ], we introduce Fermi 
coordinates (u,s) G K x [0,L] in (|3. 10.2[) (see |CdVj ) centered on the geodesic 
segment 7(0;) (i.e. u = on 7(0;)) and apply stationary phase in the it-variables 
just as in the T(H) ~ 1 case (sec [CdV, section 2]). In terms of the (it, s)-coordinates 
on M 

(3.10.3) dvol(x) — J(u, s;u>) duds, 

where, J(u, s; u)\ = O a ,p{e c{a ^ L ") [HdV] . 

In the case at hand, it is necessary to control the dependence of the phase and 
amplitude of (|3.9.6p on the parameters, lu. 

We now need the following estimates: The first is a (uniform) Hessian lower 
bound ( [CdVl Lemma 4]) which says that 

(3.10.4) VlU(u,s) >C >0. 
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The constant Co > in (|3.10.4|1 is uniform in the parameters uj € T and depends 
only the curvature pinching conditions —K\ < K < —iff. In addition, one has the 
upper bounds [Berl ICdVj : 

(3.10.5) dy (x,cjx) =O a (e CL "), d^f^u, s) = O a (e c ' L "), 

dZJ(u,s) = O a (e c " L "), 

where the estimates in (|3.10.5p are all uniform in the Fermi coordinates (u, s) S 
M. x [0, L]. In order to apply stationary phase with parameters, we need to carefully 
analyze the critical sets of the phase function fu(u, s) in (|3.10.2[) for all u> € T with 
L u <C | \nh\. We do this, by combining the estimates (|3.10.4[) and (|3.10.5[) with the 
dynamical Lemma 12.2.11 

3.11. Application of Lemma 12.2.11 As above, let j(uj) be the unique lifted 
geodesic invariant under the action of uj S T. As is well-known |B-Q| Proposition 
4.2], the displacement function fu,(x) on a negatively curved surface is strictly 
convex, except on j(uj) where it is constant: ,fu>{x) = L u when x G 7(w). The 
geodesic 7(0;) is also the critical set of fui(x), see (13.10. ip . 

Fix j = 0, [ T ^]r T ° ] Etnd u> € T and consider the corresponding summand in 
(I3.10.2p given by 
(3.11.1) 

f e^^/%(x,^)-C^^,^ t -a z r(x,z)| z= _) Xj (±r(a:,^x);r(ft)) 

xr(x, wx)' 5 77(2) dvol(x). 

We now show that, the integral in ()3. 10.2(1 is 0(H°°) unless uj — ujo, where the 
latter group element fixes the lifted geodesic 70(^0)- As before, ir : T*M — > M 
denotes the standard cotangent projection map and 7 c T*M will denote the 
bicharacteristic curve with 71(7) = 7. 

To prove the above claim, we first note that on supp (l Q (x, —duj t d z r(x, z)\ z=UJX ), 

(3.11.2) d(x,ir(j )) = d{x,To) = 0{h 5 ). 

Here, following our convention we write 70 for 70 (ujo)- Since duj^ ■ d z r{x,z)\ z=ujQX 
is the cotangent vector to 70 at x € 7r(7o), it follows from the small-scale microlo- 
calization in the ^-variables that for any a; € 70, we have that 

(3.11.3) duj 1 ■ d z r(x, z)\ z=UJX =duj t -3 + o(h d ). 

From the exponential upper bounds for the derivatives of in (|3. 10.5[) and by a 
Taylor expansion around 70, it follows that for any x € M with d(x,'fo) = 0(h s ), 

(3.11.4) duj 1 ■ d z r{x, z)\ z= ^ x = dwl ■ d z r{x, z)\ z=WoX + 0(h 5 °), 

for some < So < 5. 

On the other hand, by taking e > small enough, and again using the upper 
bounds in (|3.10.5p , it follows by an integration by parts in (|3.9.6jl in the transversal 
M-variable, that modulo 0(/z°°)-errors, one can cut off the integration to values of 
x satisfying: 

(3.11.5) d x f u (x) = 0{h s ), 
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for any < 5 < 1/2. Then using (|3. 10.51) again together with the uniform Hessian 
lower bound in (|3.10.4|) . it follows by a Taylor expansion argument that 

(3.11.6) d{ l0 {u ), 1 {u))=O{ti>i). 

Here, < <5i < <5o is yet another, possibly smaller constant. But then, in view of 
the estimate (|3.11.4|) we get that 

(3.11.7) d(7oM,7H)=0(^). 

The presence of the time cutoff \j m l|3.10.2|) ensures that there is only one j- 
summand in (|3.10.2p that contributes in a non- negligible way to k^, namely, the 
interval Ij containing L 7o . 

Finally, by possibly shrinking e > further, it follows from (|3.11.T|) and the 
dynamical Lemma \2 . 2 . 1 1 that . up to 0(h°°)-error in k^A, T), u> — luq. 

One can repeat the above argument for each j G {0, [ T ]r T ° ]} and id =^ lu g T 
separately in (|3.11.ip and sum up. By taking the exponential bounds in (|3.10.5|) into 
account and possibly further decreasing the size of e in T(K) = e\ \nh\ if necessary, it 
follows by applying stationary phase in u (see |CdVl iDonl ISunj ) that for appropriate 
a 6 > 0, 

(3.11.8) «*(A,T)= .L\ o -x(L 10 ,T(h)) ■ \ det(J - P 7o )|- 1/2 + E(h), 
where, 

f p CT(h) too \ 

(3.11.9) E{h) = 0(h a °) + O f r j = 0(h^). 

The first term on the RHS of (|3.11.9p gives the remainder produced by the sta- 
tionary phase method for lu — ujq. The second term on the RHS of (|3.11.9[) follows 
from the bound #{_7;L 7 < T} = 0{e CT{h) T(h)- 1 ), and the estimate 0{h°°) for 
each summand in (|3.10.2p corresponding to w / ui . 

Estimate for k 12 (A,T): Lemma \2 . 2 . 1 1 implies that there are no periodic geodesies 
in supp ( • (1 — C 7o ) )■ So, by repeated integration by parts in the u-variable in 
(the analogue of) (|3 . 11 . If) . it follows that 

(3.11.10) Kl2 (A, T) = O { j = <D(h°°), 
when T(h) = e\hih\. 

Estimate for k^ 2 (A,T): Here, repeat the microlocalization as in the estimate for 
Kn(A,T) near each periodic geodesic 7 ^ 70 with < L 7 < e|ln?i| separately 
and sum up. By the same argument as for «i2(A,T), all cross terms give 0(h°°) 
contributions to k(A, T). Also, we note that since the remainder term is 0(H a6 ) for 
some «6 > in (|3.11.8p . after summing over all lu 7^ ujo in nf 2 , it follows that, after 
possibly shrinking e > further, the remainder in the latter is 

0(e CT ^T(K)- 1 h ae ) = 0{h a7 ), 



for some aj > 0. This finishes the proof of Proposition 13 . 2 . ll since the 0(fi,™ 7 )-error 
is absorbed in the ©(T^-error in (jB"^ for k(A, T). □ 
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4. Proof of Theorem 11.3.11 



4.1. Application of thermodynamic formalism. The strategy of the proof of 
Theorem 1 1.3. II is to get a contradiction with Lemma l2.1.6l that gives an upper bound 
for the quantity k(X, T) denned in (|2.1.4[> . (|2.1.5|) ; that bound holds for any A, T. We 
use the formula (13.2.2j) and the estimate (|4.1.3|) below to obtain a lower bound for 
k(X, T) which holds for an infinite sequence of pairs (A, T) such that T ~ j- In In A, 
where h is the topological entropy of the geodesic flow on X. 

Consider the main term in the formula (|3.2.2p . This sum is a trigonometric 
polynomial in A, while the number of terms and the coefficients depend on T. 

We would like to choose A and T in such a way that the value of the polyno- 
mial is large. We first study the rate of growth for the sum of coefficients of this 
trigonometric polynomial as T — > oo not taking into account the oscillatory nature 
of the terms. Let 

(4.1.1) S{T)= ^ 1 



E 



V|det(/-P 7 )| 



It turns out that the asymptotic rate of growth of S(T) can be determined using 
results from the theory of thermodynamic formalism for Anosov flows: 



There exists a constant C'o > such that 



Proposition 4.1.2. 

(4.1.3) S(T) = Coe^-*^"(l + o(l)) 

as T — > co, where P is the topological pressure (jl.2.2[) and H. is the Sinai- Ruelle- 
Bowen potential (|1.2.1|) . 



It was shown in |J-P| that P {— tt) > K2/2, hence S(T) grows exponentially 
inT. 

Proof of Proposition 14.1.21 Let ^ 7 G T 7 ( )M be the tangent vector to 7. The 
Poincare map V-y preserves the unstable subspace £> and the stable subspace Ei , 
both of dimension one. The map P 7 is symplectic and has determinant equal to 
one. The eigenvalues fi, l//x oiV-y, corresponding, respectively, to the unstable and 
the stable directions, satisfy: 



(4.1.4) 



ln|/x| e ^L^KxL-y]. 



It follows from (|4.1.4|) that 



1 



det(J -V J )\ = (H - l)(l - _) = 1^1 (1 + 0(e" 



)) 



On the other hand, by definition of 7i 



so we have 



(4.1.5) 



W = exp 



jdct^-^)! =cxp 



H(G s ^)ds 



H(G s ^)ds 



(l + 0(e 



-K 2 L- 



Split the sum S"(T) into two parts: 

L y <T/2 T/2<L T <T 



Si(T) + 5 2 (T). 
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It follows from (|4.1.1|) and (|4.1.5|) that 
(4.1.6) 



S 2 (T) 



and that 



E L ~t exp 

L T GLsp,T/2<L T <T 



H{G 8 ^)ds 



(l + 0( e ~^ T / 2 )) 



H(G s ^)ds 



S 1 (T)=0\ £ 7 exp 

y L 1 ehsp,L 1 <T/2 

It follows easily from the results of [Par], [FTl (7.1)], [iVPgl p. 109] that 



S 2 (T) 



P(-H/2) 



(l + o(l)) 



and that 



5 1 (T)=o(e p (-^)?) 



This finishes the proof of the proposition 
If follows easily from Proposition 14. 1 .21 that 



□ 



(4.1.7) 



S(t) 



E 



C oe P (--)- T (l + o(l)). 



l iW -,<t Vldet(/-"P 7 )| 

Indeed, since each imprimitive geodesic is a multiple of some geodesic of at least 
twice smaller length, the contribution of imprimitive geodesies into S(T) is 

and hence could be neglected (cf. |J-P[ p. 25]). 

4.2. Preliminary estimates. We can now finish the proof of Theorem 11.3.11 As- 
sume for contradiction that Theorem 1 1 . 3 . 1 1 doesn 't hold. Then R(X) satisfies 



R(X) =0(i\n\) p ^ n ^ l -^/ h ) 



for some e > 0. Let b = P{—7i/2){l — e)/h be the exponent In A in the previous 
formula. Then by Lemma l2.1.6l we have 

(4.2.1) k(A,T) = 0((lnA) b ). 
We rewrite (14.2. lj) as 

(4.2.2) ln|re(A,T)| < 61nlnA + Ci. 

To finish the proof, it suffices to establish a contradiction with (|4.2.2[) . This will be 
done using Proposition 13 . 2 . ll and the estimate (|4.1.3[) for a suitable choice of A and 
T. In the sequel, we shall let A,T — > oo while keeping T ~ i In In A. This ensures 
that the hypothesis T < eln A of Proposition 13.2.11 is satisfied. 
Denote the main term in the asymptotics of k(A,T) by 

nx T)= y £ 7 cos(AL 7 ) x (£ 7 ,r) 

l^l^<t ?V|det(/--p 7 )| ■ 

It is a trigonometric polynomial in A. According to (|1.4.2p . without loss of generality 
we may assume that 

(4.2.3) x(L 7 ,T)> 1/2, VL 7 e(T„,T(l- £ /2)]. 
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The condition L 7 > To is not essential because it rules out only a finite number of 
closed geodesies, and their total contribution to (14.1.7P is 0(1). 

Next, we would like to choose A so that all the terms cos(AL 7 ) will be (say) 
> 1/2 for L 7 £ Lsp,To < L 7 < T. Let f(T) be the number of distinct such L 7 -s, 
and let Li,L2, ■ . . , £„(t) De the corresponding lengths. It suffices to choose A so 
that 

(4.2.4) dist(Ai J ,2 7 rZ) < 1/2, 1 < j < v{T), 
Assuming (14.2.3)) and (|4.2.4p . we get for large enough T 

ln|£(A,T)| >C 2 +ln|S(T)|. 
Using the estimate (|4.1.7I) and Proposition 13.2.11 we conclude that 

(4.2.5) ln|/c(A,T)| > P{-H/2)T{1 - e/2) - InT + C 3 . 

This formula will be used to get a contradiction with the upper bound (|4.2.2|) for 
\k(X,T)\. 

4.3. Dirichlet box principle. We next explain how to choose A so that (|4.2.4|) 
would hold. Let Mi be a large constant whose value will be specified later. Then 
by Dirichlet box principle ( [J-P| . see also |P-Rud[ iRub-Sj ) there exists 

A e [Mi, Mi2 v{T) ] 

such that (|4.2.4p holds. Hence, for any choice of M\ there exists A satisfying 

In In Mi < In In A < In In Mi + In v(T) + In In 2. 

for which (|4.2.4p holds. 

It follows from results of Margulis [MzB] that v(T) = e hT (l+o{l))/hT as T -> oo. 
Therefore, any A satisfying the previous inequality would also satisfy 

(4.3.1) lnlnMi < lnlnA < lnln(Mi) + hT - \n(hT) 
We now choose 

he 

(4.3.2) Mi = exp(exp(<xT)), where a < -. 

Then (|4.3.ip becomes 

(4.3.3) aT<lnlnA< (h + a)T - ln(hT). 

The first inequality in (|4.3.3p ensures that the hypothesis of Proposition 13 . 2 . ll is 
satisfied, implying (|3.2.2p . By the previous argument, we have shown that (|4.3.3p 
and (|4.2.5p implies existence of A such that the following estimate holds: 

(4.3.4) ln|/c(A,T)| > P(-H/2)T(1 - e/2) - lnT+ C 3 . 

To establish a contradiction with the formula (|4.2.2[) . it suffices to have 

P(-H/2)T(l-s/2)-lnT + C 3 > 61nlnA + C*i = P ^~ H '^ 1 ~ ^ ffiln A + Ci 

h 

or 

(4.3.5) lnlnA< ^ ~ £ / 2 ) T+ iL. -(C 4 -lnT) 

v ' ~ l-e P(~H/2)(l~e)- ' 

If we could show that the inequality (|4.3.3|) implies the inequality (|4.3.5p . we 
would be done. Indeed, by Dirichlet box principle there exists some A satisfying 
(|4.3.3p . and so (|4.3.5[) holds for that value of A, establishing a contradiction. The 
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linear function of T is the fastest-growing term in the right-hand side of both 
inequalities, so it suffices to compare the coefficients of T. The coefficient in (|4.3.3|) 
is equal to h + a, while that in (|4.3.5f> is equal to h(l — e/2)/(l — e). Accordingly, 
it suffices to have 

h + a < h(l-e/2)/(l-e), 
and this is ensured by the choice of a in (|4.3.2p . This establishes the desired 
contradiction and finishes the proof of Theorem 11.3.11 □ 

5. Remainder estimates in higher dimensions 

5.1. Proof of Theorem 11.5.11 In this section we assume that X is a compact 
Riemannian manifold of dimension n > 3. Recall that the Riesz mean of order k of 
a function /(A) is defined by 

(5.i.i) n k f{\) = ~£ (i - ~) f(t)dt, fc = i,2,... 



As was shown in |Safj , on any smooth compact n-dimensional Riemannian manifold 
(5.1.2) U 2 R(X) = C ai X n - 2 + 0(A"- 3 ), 

Here C is a non-zero constant depending on the dimension only, and a\ = i j„ r, 
where r is the scalar curvature of X. Note that a\ is the first heat invariant of X, 
the coefficient in the short time asymptotics of the heat trace: 

i y 1 j=0 

Therefore, if a\ does not vanish (which is the assumption of Theorem II . 5 . 1[> . 
^ 2 i?(A) » A"" 2 . Combining (j5Tl~T|) and ([5X2]) we get 

i £ \R(t)\dt >l£ (l-L^j R(t)dt = l -Tl 2 R(X) » A"- 2 . 

This completes the proof of Theorem 11.5.11 □ 
Remark 5.1.4. As follows from results of [Safj . the first Riesz mean of i?(A) satisfies 

i [ X R(t)dt = 0{X n - 2 ). 
X Jo 

At the same time, R(X) = 0(A n_1 ). Putting together these two upper bounds and 
the lower bound (|1.5.2|) one gets an idea about the amount of cancelations occurring 
when R(X) is integrated over [0, A]. 

5.2. Oscillatory error term. Following \J-P\ section 1.2] one may introduce the 
oscillatory error term R° SC (X) in Weyl's law: 

1 [^} 

(5.2.1) N(X) = - , V . flj T X n - 21 +R osc (X), 

1 ' { ' W l/2 r(|-j + i) v ; ' 

where a 3 are defined by (|5 . 1 . 3[) . The expression (|5.2.ip is not an asymptotic ex- 
pansion, however it often appears in physics literature. Such a representation is 
quite natural since it allows to separate the "mean smooth part" of the counting 
function coming from the singularity of the heat trace at zero, and the "oscillating 
part" produced by the singularities in the wave trace caused by closed geodesies. 
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We believe that using essentially the same arguments as in the proof of Theorem 
ll.3.1| one can show that the oscillatory error term on an n-dimensional compact 
negatively curved manifold satisfies: 

(5.2.2) R osc {\) = n((ln\)^ m - e ^j Ve>0. 

In order to prove (|5.2.2|) . one has to extend to dimensions n > 3 the dynamical 
part of the proof of Theorem II. 3. 11 which is easy, and to generalize Theorem ll.4.31 
which requires some work. In particular, one needs higher-dimensional analogues 
of the results of [CdVj that are used in section 3.9. We plan to carry out the details 
of this argument elsewhere. 
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